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1. Introduction 

A key test of Regge Calculus which has been studied many times in literature is the 
Friedmann universe of dust. 

The first papers on the Friedmann universe did not use simplicial complexes. The 
authors made use of non-simplicial blocks for which they had to give some other 
information besides the edge lengths. This information was deduced by symmetry. 
Collins and Williams (1973) used three tridimensional regular polytopes to describe 
the spatial sections with topology 5"^ of the Friedmann universe of dust. Then Brewin 
(1987) wrote in the discrete formalism of Regge Calculus the expression of the action. 
Following Barrett et al (1997), for a regular 3-polytope having n vertices, matter is 
placed on the timelike edges joining a vertex i of the polytope with its evolute i' of 
the next spatial section, finding: 

b i 

where the first sum is over the bones b of the manifold and Eb is the defect of the bone 
b with area Ab. The second sum is done over every vertex i of the first spatial section. 
Since a regular polytope having n = 120 has been chosen as spatial section, a mass 
M/120 is placed in each vertex i. Besides the regular polytopes, Brewin also used 
other non regular tridimensional polytopes. The evolution of all his polytopes, using 
Regge equations, stopped before the volume of spatial 3-section became zero: a point 
was reached where the equations had no solution. He saw that this point approaches 
the limit value (spatial section of null radius) as the polytope approaches the 3-sphere 
of the continuum solution. To overcome the obstacle, Brewin continued the evolution 
reverting to differential equations. 

The first authors who used a simplicial complex for this metric were Barrett et 
al (1997). They evolved the Friedmann universe of dust using the Regge equations 
found by extremizing the action (1): 

EdAijk n 
k ^^ 

where the sum is done over all vertices k joined to the edge [ij] of length Uj. The bone 
[ijk] has Eijk as defect and Aij^ as area. They made use of the evolutive scheme found 
by Sorkin (1975), but in our opinion their simplicial complex is not built correctly 
(see subsubsection 3.2.1). However they too found the point of stop and saw that 
its position is independent of the timelike interval between two consecutive spatial 
sections. The authors advanced the idea that the evolution should be continued by 
making spacelike the edges of the lattice which were timelike before the stop. However 
they did not pursue the idea. 

In the present paper we propose to solve the problem of the stop of the evolution 
and to understand the meaning of the ending point. The paper is organized as follows. 
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In section 2 we embed the Robertson- Walker metric in a five-dimensional Lorentzian 
space. We will find it useful in order to understand the problems related to the 
evolution of our lattice. Then we describe our lattice and carry out the first evolution 
in order to create the initial simplicial sandwich (section 3). Section 4 contains 
the evolution of the initial sandwich and the results obtained. The description of a 
correct evolution algorithm for our simplicial complex based on Regge equations and 
the Sorkin evolutive scheme is also presented here. The nature of the singularities 
belonging to a general closed Robertson- Walker metric is discussed in section 5. 
In section 6 we propose a deeper reason why our simplicial approximation of the 
Friedmann universe ceases to exist before the volume of the spatial 3-section vanishes. 
We find that the stop condition is caused by a novel singularity of the metric. We 
fit the simplicial solution with a continuous metric and study its stress-energy tensor 
(section 7). Finally, our conclusions are reported in section 8. 

2. Embedding into a 5-dimensional space 

We shall find convenient to embed a space-time with closed Robertson- Walker metric 
into a 5-dimensional Lorentz space. This will be of help for understanding the 
behaviour of the Regge equations when applied to the evolution of a 600-cell space 
section. 

The metric in consideration can be written in this way: 

dx" dx^ = -dt^ + R^ [dx^ + sin^ X {d^^ + sin^ ^ d^^) ] , (2) 

where R depends only on the universal time t. 

The variety with metric (2) can be embedded into a 5-dimensional space endowed 
with the following metric: 

da;" dxa = -d^;^ + dw"^ + dx"^ + dy'^ + dz^ 

by putting 

w = R cos X 

X = R sin X sin t? cos </? 

(3) 
y = R sin X sin 1) sin </? 

z = R sin X cos t?. 
Equations (3) imply 

dx" dxa = -dv^ + dR^ + R^ [dx^ + sin^ x (d-i?^ + sin^ ^ dv?^)] . (4) 

In order that metrics (2) and (4) may coincide it is necessary that 

df = dv^ - dR\ (5) 
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If we impose that R is a, function only of the time variable v then we have 

dR= ^dv, 
av 

so 

dt' = dv' - (^)'d^' = (1 - R') d^', 

from which it follows that R^ < 1. In this case 

dt = (1 - R^f^ dv. (6) 

So the metric becomes: 

dx" dxa = - (1 - R^) dv^ + R^ [dx^ + sin^ x (d^?^ + sin^ {} df/:^)] . (7) 

In the following we will refer to v as "outer time." We will be mainly interested in 
the case where R{v) is an even function, so that t; = is an instant of time symmetry. 

3. The 600-ceII and the construction of the initial sandwich 

3.1. The spatial section 

We will approximate the hypersphere with the last tridimensional regular polytope, 
the 600-cell. It is a simplicial complex constituted by 120 vertices and 600 tetrahedra. 
The Priedmann metric has an instant of time symmetry where we can put a spatial 
section with a null extrinsic curvature. In Regge Calculus, for a spatial section whose 
extrinsic curvature vanishes, the following relation must hold: 

£io- liA = 167r 



^-^'^'^^ 120' 



J 



where the sum is done over all vertices j joined to vertex i by an edge [ij] of length 
lij and defect Sij. There are 5 regular tetrahedra around a single edge, and the defect 
of each edge is: 

£ = 27r — 5 arccos(l). 

Since every vertex is joined to 12 vertices which make an icosahedron, and the 600-cell 
is regular, calling /q its edge, we have 

90 

M^—elo. (8) 

IT 

In order to have a criterion to compare the results of Regge Calculus with those of 
General Relativity, we choose that the 600-cell should approximate a 3-sphere having 



The Friedmann universe by Regge Calculus 5 

the same volume. Thus the radius Rq of the 3-sphere, in the instant of time symmetry, 
will obey the relation 

73 

27r2 Rf, = 600 " 



•-o 



6^2 



or 



k — CRo 



where 



7rV2 
50 



1/3 



3.2. Creation of the initial sandwich 

We call "initial sandwich" the space-time between the initial spacelike section (Eq) 
and the next one (Si). 

3.2.1. The classes of vertices: a, /?, 7, 5, £. Tuckey (1993) and Barrett et al (1997) 
say the 120 vertices of a 600-cell can be subdivided into 4 classes. In each class there 
should be 30 vertices not joined to one another by an edge. So all vertices of a class 
could be evolved simultaneously. 

But this cannot be true: it can be shown as follows. Let us consider a vertex V: 
we know it is surrounded by 12 vertices Wi making an icosahedron. These vertices 
cannot be of the same class of V because every vertex Wi is joined to V by an edge. 
But in an icosahedron it is not possible to choose more than 3 vertices not joined to 
one another, so that these 12 vertices must belong to at least 4 distinct classes. 

To overcome this obstacle let us reason as follows. Let us consider V and let it be 
a vertex of class a. Let us try to subdivide the vertices Wi into 4 groups of 3 vertices 
not joined to one another. We will name these groups with the letters /3, 7, 5, e. This 
can be done by looking at figure 1. 

Let us try to subdivide the vertices of the 600-cell into 5 classes of 24 vertices 
each. Every class must consist of vertices not joined together. Using the classification 
of the vertices of the 600-cell given by Coxeter (1973), we find that 

Cl = 1^07 ^107 ^207 ^307 ^407 ^507 Bj, Bij, B27, -B37, -B47, -B57, 

Ci, Oil, C21, C31, C41, C51, Ds, D18, D2S, -D387 -D487 D^s}- 

The other classes are obtained by adding 2, 4, 6, 8 to the subscripts. It is also 
possible to show that the 24 vertices of each class make a 24-cell regular polytope 
inscribed in the 600-cell. Furthermore we can see that every vertex is actually 
surrounded by 12 vertices equally subdivided among the other classes. So an a- vertex 
is joined to 3 /9-vertices, 3 7-vertices, 3 5-vertices and 3 e-vertices. It should be noted 
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that the 3 vertices of the same class joined to a same vertex are not equivalent. For 
example the 3 /9-vertices around the vertex a in figure 1 have different neighbouring 
vertices. In fact they are surrounded by three sets (of 5 vertices each), which are not 
equal to one another. By the edge [a(3o] we will call the edge [a(3] surrounded by 
the set of the following vertices {7, 7, 5, 5, e}. The edge [a(3i] is surrounded by the set 
{7,7, 5, £,£} and [a/32] by {7, 5, 5, £, e}. As we shall see later, this inequality among 
the three edges [af3i] (z = 0, . . . , 2) will give birth to three different Regge equations 
when an a- vertex is going to be evolved. 

Because of all these facts the evolution of a 600-cell consists of five steps: the 
evolution of each class of vertices (a, /?, 7, 5, e). 

3.2.2. The initial sandwich and the instant of time symmetry . Equation (8) has been 
found by choosing a spatial section having a null extrinsic curvature. In the continuum 
this means that the orthogonal timelike geodesies, which go out from this section, are 
locally parallel. In a lattice too, this characteristic should remain. So we can say 
that the spatial section, found by putting /q equal to the value given by equation (8), 
cannot be a section of time symmetry. In fact if this were true, the edge lengths in 
the next section would be smaller: V = li < Iq. Then the timelike edges which are the 
geodesies of the masses could not be parallel. 

To understand better what happens in this situation let us embed the 600-cell into 
the 5-dimensional lorentzian space described in section 2. Since the 600-cell can be 
inscribed in a 3-sphere, we can look at its evolution as a sequence of 3-spheres each 
having a radius R{v), where v is the outer time of the 5-space. Then we put 

a = {Vj Rua) 

where Ua is a unit spacelike 4-vector, 

Ua = (cOSXa, siu Xq, siu Z?^ COS (/7q,, siu Xq, siu 1?q siu (/7q , sin Xa COS 'da) J 

pointing to the direction of a. There will be one such vector u for each vertex of all 
600-cell spatial sections. Let us call Eq the section corresponding to the instant of 
time symmetry v = 0, and Ei the next one, at time vi, then 

a = (0, RoUa) (3 = (0, RqU/s) 

a' = (vi.RiUa) (3' = (vi.RiUp). 

So we have that [aa'] is not orthogonal to [a/?] , unlike what happens in the continuum. 
In fact the straight lines 

X = a + p [aa'] 

^ = (3 + q[(3(3'] 

are never parallel and cross forp = q = Rq/{Rq — Ri). 
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Nevertheless we can make an attempt in order that our lattice may have a 
behaviour similar to that of the continuum. Let us now put Eq and Ei symmetric 
about the instant of time symmetry v = 0, then vq = —Vi. Now we have 



a = {-Vi^RiUa) 


13 = {-vi^Riu^) 


a' = { Vi,RiUa) 


P' = { v^^R^Uf,). 


Thus we find that 




a'P'] = [aP], 




aa'] = IPP'] = (2 


z;i,0,0,0,0), 


and [aa'] is orthogonal to [a/3 . 


Therefore we obtain 


a' (3] ■ [a' (3] = [aa'[ 


■ [aa'] + [a(3] ■ [a(3 = 



.4 ,,2 1/2 ^ _o-2 , /2 

'I C]^ -t- (,Q ~ /q -t- tg, 

where Tq is the interval of proper time between Eq and Ei. In fact we have also tq ~ Av 
because R = for v = 0. 

Now all timelike edges (geodesies of the mass points) are orthogonal to the spatial 
section and parallel to one another, as we would expect if the extrinsic curvature 
vanished. Equation (8) does not give the edge of the time symmetric spatial section, 
but the edge of two consecutive 600-cells equal and symmetric about v = 0. All other 
spatial sections will also be pairwise symmetric about v = 0. 

So we know the squared lengths ' of all edges of the initial sandwich and are able 
to build it at once without using any Regge equation. 

3.2.3. The initial sandwich using Regge equations. To test our assumptions we evolved 
the section Eq making use of the Regge equations for the edges lying between Eq and 
El, according to the Sorkin evolutive scheme. In order to allow an easier comparison 
between our results and those of Barrett et al (1997) we have chosen a lattice having 

M = 10.202 . . . (in arbitrary units), 
/o = 2.774..., 

To = 0.0102. 

In fact the authors considered the Friedmann universe having maximum radius Ro 
and mass M^, defined by 

M, = fg{t) dV = 2^ qR' = ^ Ro. 

Choosing M* = 10 they found 

lo = CRo = ^M, = 2.774... 

' We call "squared length" of an edge [ij] the real quantity s^ — [ij] ■ [ij], whereas its length is the 
non- negative number l^j — \s^A^'^ 
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In the end they found, using equation (8), the value 10.202 ... for M. 

We have aheady remarked that the evolution scheme does not preserve the exact 
symmetry of the initial 600-cell. We expect however that, if the evolution is properly 
done, the symmetry is approximately respected. In order to test this, we examined 
those variables which are bound to remain equal if the symmetry holds and computed 
their averages and standard deviations. These variables are: the squared length d'^ of 
the oblique edge joining Eq and Ei and the squared length /' of the edge of the new 
600-cell forming the Si section. 
We obtained 

ll = 7.693 799 901 383 04 
/2 _ T-2 = 7 693 695 861 383 04 

d^ = 7.693 695 861 383 01 a^^ = 1-348 • 10"^^ 
f = 7.693 799 901 382 97 a^,^ = 5.767 • 10"^^ 

So, as we expected, the same value /q occurs on the first two consecutive surfaces. 
Furthermore we see that the relation (i^ = /q — Tq holds up to the 14th significant 
figure. 

As we already remarked, the other surfaces, to be obtained later, are symmetric 
about these two. This last result has been found by Barrett et al (1997) too. The 
authors also found a second time-symmetric solution, in which the maximum spatial 
length is reached on a single hypersurface and the other hypersurfaces are symmetric 
in pairs around it. However this new maximum is slightly larger than that given by 
equation (8). These are all acceptable models and are indistinguishable, on the scale 
of time we used, from ours. 



4. The evolution of the initial sandwich by Sorkin evolutive scheme 

Given the squared lengths of all the edges lying between and on two spatial sections 
Efe_i and E^, in order to create a new simplicial sandwich, i.e. the space-time between 
Efe and Efc+i, we will follow five steps. 

1. To evolve the vertices of class a' of E^. We put a timelike edge (of negative 
squared length Sa"a' = —T~k) joining a' to its evolute a" of E^+i. We link a" to 
the 12 neighbouring vertices of q;"s entourage in E^ through spacelike edges (to 
ensure this, |sa"Q'| should be small enough with respect to the other edges). 
There is one Regge equation for each "internal" edge. Then we find Regge 
equations around [q;'/3j'], [ct'7i], [ct'^i], W^'i\ (i = . . . 2) and around [q;"q;']. We 
have as many equations as unknowns. The problem looks completely determined. 
Things are not so easy, however. 

As widely shown in literature, on one side we must save the "lapse and shift" 
freedom; on the other side we expect Regge equations to be not independent 
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because of the contracted Bianchi identities. Therefore we must ignore 4 
equations and fix 4 of the unknowns by proper lapse and shift conditions. As 
condition of lapse we will choose at our convenience 

_ 2 

Sa"a' — T"fc 7 

and as conditions of shift the three constraints 

We then cancel the 4 equations for [a" a'] , [a'(3[] , [a'f]!^] , [ct'7i] • 

As already remarked in subsubsection 3.2.1, equations over edges of the same 

class, e.g. [a'S'j] are different. In fact it could be seen that [a'^g] is linked to 8 

vertices (one d, one £, two /3', two 7', one e' and one a") whereas [a'd[] is linked 

to 9 vertices (one 5, two £, two /3', one 7', two e' and one a"). 

There are still the following 9 unknowns: 

-Sa"/3^, ■Sa"7^, ■Sa"7^ , Sa"S'-, Sa"e'-- 

Nevertheless, we can say that the 12 slanting edges belong to 4 different classes, 
corresponding to the class of the second vertex. It is natural to require, in order 
to save the symmetry, that edges belonging to the same class be of equal length. 
So the unknowns reduce to 4: 

Sa"/3', Sa"Y^ Sa"5', Sa"s'j 

after putting equal all connections [a"d'] , all connections [a"e'] and Sa"-y' = Sa"-y!^ ■ 
So the equations must reduce to four, for instance: 

[a/3o], [a'^o]. ["'To]. ["'4]- 

In the following we shall verify that the other equations we have not used are 
satisfied anyhow. 

To evolve the vertices of class /?'. We join (3' to (3" and link /3" with the 3 vertices 
of class a" of its entourage in E^+i and with the 9 vertices (of classes 7', d' and 
e') of /9"s entourage in E^. By symmetry, we may reduce the new edge lengths 
to the following: 

Sl3"Y, Spi'S', Spile', Sa"/3", 

since we put as lapse s 1311/31 = — r|. The equations are 

[(3'^l [(3'dl 1(3' el [a"n 
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neglecting the equation for [(3" (3'] because of Bianchi identities. 

3. To evolve the vertices of class 7'. We join 7' to 7" and link 7" with the 6 vertices 
(of classes a" and (3") of its entourage in E^+i and with the 6 vertices (of classes 
5' and e') of 7"s entourage in E^. Now we have five kinds of internal edges: [7"7'] , 
[q;"7'], [/3"7'], [7'^'] and [y'e']. These give rise to five distinct Regge equations. 
But because of the Bianchi identity only four of these can be used to find four 
lengths (neglecting the one for the timelike edge), whereas syry must be freely 
given, like Sa"a' and put equal to — r|. Therefore we take as unknowns 

Sys', Syi^i, Sa"y, Sp'ry. 

4. To evolve the vertices of class 6'. We join 6' to 6" and link 6" with the 9 vertices 
(of classes a", (3" and 7") of its entourage in S^+i and with the 3 vertices of class 
e' of J"s entourage in E^. By symmetry, we impose that the new edge lengths 
are: 

S5"e', Sa"S", Sf3"S", Sys", Ss"5'- 

We have only five distinct Regge equations, around [S"d'], [a"d'], [f3"d'], [7"^'] 
and [S'e'] . One of these, that for the timelike edge, will not be utilized because 
of the Bianchi identity. So the remaining four can be used to find four lengths, 
whereas ss"S' is given, like Sa"a' and equal to — r|. 

5. To evolve the vertices of class e'. We join e' to e" and link e" with the 12 vertices 
(of classes a", (3", 7" and 6") of its entourage in "Ek+i- Now we have five kinds 
of internal edges: [e"e'], [a"e'], [/3"e'], [y"e'] and [S"e']. These give rise to five 
distinct Regge equations. But because of the Bianchi identity only four of these 
can be used to find four lengths (neglecting that for the timelike edge) , whereas 
Se"e' must be freely given, like Sa"a' and put equal to — r^. Therefore we take as 
unknowns 

Sa"e"j Sp"e"j Ss"e", Sy's"- 

It should be noted that every evolution consists in solving five systems of 4 
equations in 4 unknowns each. Furthermore the simplicial sandwich does not have 
the same grade of symmetry as the initial spatial section, where all vertices were 
equivalent. The evolutive algorithm has broken the symmetry among vertices a, (3, 7, 
5, e. 
The evolution of E^ is now finished: the new sandwich has been built up. 

4.1. Results of the evolution 

We evolved the initial 600-cell by imposing that the edges of the same type, like [a'/9j'] , 
were be all equal. In this way only four unknowns and four equations remained when 
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each vertex was evolved. Actually the different variables should be only two, since 
symmetry requires that all edges [1;"^'] and [t'"zy"] are equal as well. This was done by 
Barrett et al (1997). By allowing the edges joining vertices belonging to two different 
classes to be possibly different one from another, we introduce new degrees of freedom 
in the lattice. We have verified that this enhancement of freedom causes no problem 
in the evolution. 

4.1.1. Behaviour of Rk as a function of the universal time t. Using the evolutive 
algorithm we evaluated, for each section, the edge (or better the mean edge) of the 
600-cell, Ik- Then we calculated the quantity 

„ Ik 

ID "■ 

where Rk represents the radius of the 3-sphere equivalent to the 600-cell of the k-th 
spatial section (E^) generated by the k-th evolution. We have that Ik = Ikit)^ where 

fc-i 

Tj is the interval of proper time between Ej and Ej_|_i, and tq = 0.0102. As better 
explained in the appendix, we did not choose equal proper time intervals. Instead we 
took Ti = pli, with p ^ 3.68 ■ 10"^. 

The results of this evolution can be seen in figure 2. 
Figure 2 shows that the evolution stops at a finite value of the radius (i?m ~ 1.2) 
which is approximately \ of the initial value {Rq = 4.244...). For R < R^ Regge 
equations give no solution. 

As it can be seen in the appendix our evolution steps were done with a constant 
interval of conformal time: equation (Al) says that Ar] ~ 2.4 ■ 10~^. Actually the 
critical point occurs after k^ ~ 750 evolutions, or rym ^ 1.8. But, since rj should reach 
TT, we should have done a number 



' 1 



Ar] 



k^^ = ^^ 1300 



of evolutions. This behaviour of our "platonic" polytope was also found by Barrett et 
al (1997) and Brewin (1987) who also studied not-simplicial complexes that, compared 
to the 600-cell, better approximated the 3-sphere. For these polyhedra the graph of 
R{t) is nearer to the cycloidal one of the Friedmann metric and the point of stop is 
placed in ?7 > ?7jn. 

4.1.2. Behaviour of the two a. During the evolution, the two a ((t^2, 0-^,2) remain 
below 10^^^. In proximity of the critical point these variables grow up quickly. This 
means that the lattice is twisting more and more. 
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4.1.3. Bianchi identity. We calculated, for each evolution, the residual fa"ai of the 
Regge equation for the timelike edge [«"«'], 

Y^ dAaiia'k "^ 7,# 

Ja"a' — / , £a"a'fc —;^, TT ^^^ • 

^-^ Ola"a' 15 

k 

This equation had not been utilised to build up the spatial section E^+i because of the 
simplicial version of the contracted Bianchi identity, i.e. the energy conservation in the 
form of the Kirchhoff-like law. The smaller fa" a', the better this identity is respected. 
It should be noted that the Kirchhoff-like law always exactly holds because the current 
of matter along the timelike edges remains the same after each evolution. Instead, 
Regge equations represent the Einstein tensor only approximately (Barrett (1986)). 
Our results can be seen in figure 3: we see that fa"a' is going to "explode" in proximity 
of the critical point. 

The other equations, which are not utilized during the evolution, are all over 
slanting and orizontal edges. We found that they are satisfied with greater accuracy, 
since their residuals are all less than 2 ■ 10^^^. 



5. Analysis of the critical point 

What is the nature of the critical point? Brewin (1987) gives a partial answer by 
noting that the following relation occurs: 

A/ ^ 

r 

i.e. the collapse becomes so fast that the vertical edges are to become spacelike. 
Unfortunately in our model the above fraction reaches a not easily interpretable value 
(approximately 2.6). So what is the true meaning of this condition, and how can we 
continue the iterations (if possible)? 

To answer these questions, it is useful to study more deeply the characteristics 
of the homogeneous and isotropic metric with topology (and geometry) 5"^. The 5- 
dimensional embedding described in section 2 will be of use here. 

5.1. Singularities 

In order to study the singularities of metric (7) we can compute the quadratic Riemann 
invariant: 



12 

Q = R"^ ^ RafS^S 



R^l-R^y 



R^ J_ 

(l-i?2)2 R\ 



(9) 



In fact a necessary and sufficient condition for a singularity is that Q becomes infinite. 
A glance to equation (9) shows that the singular points are those values of v where 
\R\ = 1 or R = or \R\ =+00. Let us discuss the possible cases. 
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1. The points R = 0. These coincide with the infinite contraction of the universe 
when the volume of the spatial section goes to 0. In these points, as will be shown 
later, density reaches infinite values. By a volume-vanishing (VV) singularity we 
will mean one of these points. 

2. The points R = ±1. If we are not in the first case we can say that the radius of 
the space section increases or decreases with speed tending to that of light. If one 
tries to extend the manifold beyond these points, and if |-R| increases further, the 
metric becomes positively defined. So it can no longer represent a manifold of 
General Relativity. This is tantamount to saying that the spacelike sections are 
not causally connected: the points of two different spatial sections can be joined 
only by spacelike geodesies. So these points will be said to be causality-breaking 
(CB) singularities. Actually, in a rigorous sense a CB singularity is a border of 
the manifold, and the extension is meaningless. 

3. The points R = ±oo. If we are not in the previous two cases the singular 
behaviour is shown by the geodesies deviation, which becomes infinite. We will 
refer to these points as geodesies-deviation (GD) singularities. 

These singularities can also occur together in the same point: an example is the 
hypothetical universe whose radius has the form R{v) = v — A v^^"^ (the critical point 
t) = is a VV-CB-GD singularity). 

5.2. Study of the stress-energy tensor 

From Einstein equations for the closed Robertson-Walker metric we have for the 
diagonal components of T^^ in an orthonormal basis: 

Stt 1 fdR\^ 1 



3 ^ R^\dt J ^ R^ 

1 d^R 1 fdR\^ 1 



R dt2 2R^ \dt J 2i?2- 
By expressing the time derivatives in terms of the outer time we find 

3 1 

Q 



P 



^^R^l-R^) 

(10) 
1 2RR+1-R^ 



Stt R^{l-R^y 



5.2.1. Range of validity of the stress- energy tensor. It is well known that the stress- 
energy tensor must satisfy 

^> \p\- 
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In our case 

-2<-?^<l. (11) 

l-i?2 

In general it can happen that condition (11) be not satisfied for values of v for which 
the metric still has a "geometrical sense." In fact a CB or GD or CB-GD singularity 
does not satisfy equation (11). Instead a VV point presents no problem. It should 
be noted that we cannot become aware of these new "physical" matter-limiting (ML) 
points, wherein g = \p\, by looking only at the expression of the metric. 

5.3. Friedmann metric 

Since dt — Rdt], from equation (5) we obtain 



V 



f[R^-{dR/driYf^'dv. 



For the Friedmann metric we have: 

V — ZKo sm — 
and 

2 

R = Rocos^- = Ro- -^. 
" 2 4Ro 

Metric (7) takes the following form: 

dx" dx^^-^ dv^ + R^ [dx^ + sin^ x (d^^ + sin^ ^ d^'^)] . 

Rq 

Calling Vc the value of v where i? = — 1 and v^^ the positive value of v where 
i? = 0, for the Friedmann metric these two singularities occur at the same time, 
i.e. v^ = Vc = 2i?o- So, in the 5-space, in correspondence of the big crunch 
the speed of contraction equals the speed of light, i.e. the big crunch is a VV-CB 
point. This behaviour is shared by all matter satisfying p = w g, where to is a 
constant, — | < w < 1, since then R^ = 1 — A^ R^'^^'^, A a constant. Seeing that 
R — — 1^2(1 + 3w)R^'^, in Zel'dovich's interval, < w < 1, no GD singularity is 
present beside the other two. 

6. The problem of the 600-cell 

Let us now reconsider equation (5). For each evolutive step we know At = t and 
AR = A//C. Then we can find Av by calculating Av ^ (At^ + AR^f^^. So we 
determined i? as a function of v (figure 4). 
We have also made a picture of AR/ Av ~ dR/dv as a function of R in figure 5. 

We can see that this derivative in the last iterations tends, in absolute value, to 
1. In this sense the 600-cell meets a CB singularity before its volume vanishes. So we 
expect that for the 600-cell Vc < v^. This leads us to believe that there is no point 
in trying to continue the evolution beyond this point, perhaps by making all edges 
spacelike as suggested by Barrett et al (1997). 
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7. Study of a "new" class of solutions of Einstein equations 

We have just seen that the behaviour of the evolution for the 600-cell looks somewhat 
different from the Friedmann solution, in its reaching a CB singularity before the 
vanishing of R. So it appears expedient to study a generalisation of Friedmann metric, 
also having that property. For the 600-cell an empirical fit of R{v) shows that a good 
fit is given by 

R{v) = i?o - ^, with a^ ^ 1.128, 

4-Ro 

whereas a = 1 is the Friedmann metric. 
It is useful to re-define 

def 2Rq 

Vc = -^, 



def 2ixo 

a 



Vn = = (^^r. 



so 



Vc 
Vr. 



• For < a < 1 we have that —Vq^ < v < v^ and the singularity is of VV-type. 

• For a > 1 we have that —Vc < v < Vc and the singularity is of CB-type. 

• For a = 1 the solution is the Friedmann metric. So this case is a watershed 
between two different behaviours of the general metric. 

It is not difficult to find the expression of t as a function of v. 

1 /9 

t = \vc arcsin(i;/t;c) + | f [l — {v/vcf''\ 

and to show that the graph of R{t) is still a cycloid, like for Friedmann metric, but 
scaled in i? by a factor l/a? and translated upwards by \vc (a^ — 1). 

7.1. Behaviour of the metric 

Condition (11) implies that the following inequality must be satisfied: 

2^-2'^'^^ 2 1 2 12 1 I 

V <v = v^ — ^v^\a —1\. 
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If a < 1 this condition does not carry any problem because v'^ > v^. If a > 1 the ML 
point V is the first one the metric meets in its evolution. In fact 

~2 2 12/2 i\^2 

Note that in order to have -D^ > the value of a^ must not exceed 4. So the conditions 
for the stress-energy tensor reduce the interval of definition of a to < a < 2. If 
a < 1, the metric ends in a VV singularity reached the more softly the smaller the 
value of a. If a > 1, the metric reaches a ML point before the CB one. It should be 
also pointed out that the Friedmann metric (a = 1) is more "pathological" because 
its VV singularity is also a CB point, since v^ = Vc- 

For this metric it is also possible in principle to write down an equation of state 
of matter, but of course we do not attach a physical meaning to this result. 

8. Conclusions 

We have seen that the evolution of the 600-cell does not describe the Friedmann 
universe well. Instead we can think of it as the evolution of a different type of matter. 
But where does such a difference come from? In our opinion, one should consider the 
following two points as playing an important role in answering this question. 

First, the 600-cell is a good approximation to a 3-sphere, but it is still too 
"rough." For every approximating method, Regge Calculus included, the smaller 
the step interval the better the fit. When this interval reduces, the two solutions 
(numerical and analytical) tend to coincide. In Regge Calculus, in order to get a 
better approximation, one should increase the number of tetrahedra of each space 
section, i.e. the 600-cell should be substituted by another non-regular polytope nearer 
to a 3-sphere (Brewin (1987)). 

The second aspect is deeper since it is closely related to the grounds of the 
Regge Calculus itself. We have seen that Regge equations give rise to a different 
evolution of the universe, ending in a CB point. This happens because they are only 
an approximation of Einstein equations. The point has been discussed by Brewin and 
Gentle. It can be said that Regge equations are a good approximation to differential 
equations rather different from Einstein's. It is wholly possible that such equations 
have solutions which behave in a qualitatively different way from those we are looking 
for: for instance, the approximate solution may exhibit singularities not belonging 
to the correct solution. Nevertheless, it remains true that the smaller the truncation 
error, the closer the behaviour of the numerically approximated equation to that of 
the exact solution. 

Thus an idea comes to mind: could we modify Regge equations in order to have 
a smaller truncation error? It has been shown (see Brewin and Gentle, Gentle and 
Miller (1998)) that the solutions of Regge Calculus are, in general, expected to be 
second order accurate approximations to the corresponding continuum solutions, i.e. 
the truncation error for lattices obtained by Regge Calculus is O (/'^) . Thus a research 
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program could be: to find a modification of Regge equations giving a truncation error 
0{l^) or better. 

Appendix. Numerical topics 

We did not use a fixed r^ for all evolution steps. Instead we imposed a constant ratio 
between Tk and Ik- This in order to ensure that the Courant criterion continues to 
be respected also when /^ becomes very small. We remember that the 600-cell should 
approximate a Friedmann universe having the same volume. So the radius Rk of the 
k-th 3-sphere is linked to the edge Ik of the 600-cell of the spatial section E^ of the 
k-th evolutive step by the relation Ik = ( Rk (C defined in subsection 3.1). Then 

^ = p'^^^ 3.68 ■ 10-^ 
h lo 

This implies that 

^^"^ = ^' = ^'^"2.4.10-^ (Al) 

K Kk 

To solve Regge equations, a system of non linear equations, we use the Newton- 
Raphson method. To begin the iterations we have to give values to the unknowns in 
order to start the iterations. We call them "trial values" and they are marked by a 
tilde. 

To give a trial value to the variables ^'''^Sa"v' (the squared lengths of the oblique 
edges that join a" to the vertices v' ^ ol of E^) we can use the following argument. 
Keeping the notations of subsubsection 3.2.2 on, we can write 



bra" 1)' 



a' a"] ■ [a' a"] + [a'/?'] + [a'/?'] - 2 [a' a"] ■ [«'/?'] 
-rl + ll~2[a'a"]-[a'n 



where 



[a' a"] ■ [a' 13'] = (Ai^^, {Rk+i - Rk) u^) ■ (O, Rk {up - u^)) = 

d R 

= Rk {Rk+l — Rk){Ua ■ W/3 — 1) ~ ~~a7^^^^ TkRkCj 

being Tk the interval of universal time between the sections E^ and E^+i, Rk the radius 
of the 3-sphere circumscribed to the 600-cell and 

^ def ^ ^ ^ 3 - V^ 
4 = l-Ua-U/3= , 

found by calculating the angle between two vectors pointing to any two neighbouring 
vertices of the 600-cell (Coxeter (1973)). Furthermore we can see that Ik = 4>Rk^ 
where (f) is the golden ratio, \ {\fb — 1). So we are able to write 

W5^„^, = _^2 j^i2^_2i^rklk {k/lk - if^. 
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Once the vertices a' have been evolved we put the quantity ^'^^Sa"v' as trial value 
for all other oblique edges, ^'^^sp'^v/, ^''^s-yi'v', ^''^ss"e'- At the moment of the evolution of 
the vertex (3' we have to give another trial value to the unknown '•^-'sq/z/j//. 

In the Friedmann metric we would have 

R{r])=Rocos\r]/2) 

so, since / = C -R, expanding to the second order in r, we can calculate an approximate 
value for the edge of the 600-cell of the section E^t+i. In fact we have 

Ik+i ^Ik- n C {lo/h - 1) - i C^ ^0 (Tk/lk) ■ 

After the evolution of /^'-vertices the quantity ^'^^Sa'/p" is determined and we use it to 
give the trial values for all variables representing the edge of the 600-cell of the section 
Sfc+i, i.e. ' 's^'/y'/, ' 'ss"v", ^ 'S£"v"- 
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Figure 1. Vertices around an a- vertex 
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Figure 2. Evolution of tlic radius i? as a function of the universal time t according to 
Regge Calculus and General Relativity 
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Figure 3. Residual f^ir^r vs. t 
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Figure 4. Behaviour of i? as a function of v. 
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Figure 5. Behaviour of |A_R/Aw| ss \R\ as a function of v for tlie 600-cell 



